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within the dynamical mean-field theory and using the geometrically averaged ("typical") local 
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with each other stabilizing the metallic phase against occurring one of the insulators. The Mott and 
Anderson insulators are found to be continuously connected. 
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I. INTRODUCTION 

Motion of quantum particles can be suppressed or even 
destroyed by Coulomb interactions and disorder, which 
are the driving forces behind a metal-insulator transition 
(MIT). The Mott-Hubbard MIT is caused by Coulomb 
correlations in the pure system, i.e. without disorder^ 
The Anderson MIT, also refereed to the Anderson local- 
ization, is due to coherent backscattering from randomly 
distributed impurities in a system without interaction.-^ 
The properties of real materials are strongly influenced 
by both interaction and randomnessi^ It is therefore a 
challenge to investigate quantum models where both cor- 
relations and disorder are simultaneously present 

The Mott-Hubbard MIT is characterized by opening 
a gap in the density of states at the Fermi leveli^ At 
the Anderson localization the character of the spectrum 
at the Fermi level changes from continuous one to dense 
pure point onei^ It is plausible that both MITs could be 
detected by knowing a single quantity, namely, a local 
density of states (LDOS). Although the LDOS is not an 
order parameter associated with a symmetry breaking 
phase transition;^ it discriminates between a metal and 
an insulator, which is driven by correlations and disorder. 

In a disordered system the LDOS depends on particu- 
lar realization of the disorder. Then the entire probability 
distribution function of the LDOS is required to know,^^ 
which is a very demanding task. Instead one could use 
certain moments of the LDOS. This however is insuf- 
ficient because the arithmetically averaged LDOS (first 
moment) stays finite at the Anderson MIT.^^ It was al- 
ready pointed out by Anderson* that the "typical" values 
of random quantities, which are mathematically given by 
the most probable values of the probability distribution 
functions, should be used to describe localization. The 
geometric meauMi^ gives an approximation of the most 
probable ("typical") value of the LDOS and vanishes at 
a critical strength of the disorder, hence providing an 
explicit criterion for Anderson localizationiSiiSiiSiii 

Theoretical descriptions of the MIT has to be non- 
perturbative if no long-range order exists on either side 
of the transition. In such a case there is no obvious or- 
der parameter and no Landau type functional available. 
A non-perturbative framework for investigation of the 



Mott-Hubbard MIT in lattice electrons with a local in- 
teraction and disorder is given by dynamical mean-field 
theory (DMFT)iS*i2i2a*2i If in this approach the effect 
of local disorder is taken into account through the arith- 
metic mean of the LDOS^ one obtains, in the absence of 
interactions, the well known coherent potential approxi- 
mation (CPA),^^ which does not describe the physics of 
Anderson localization. To overcome this deficiency Do- 
brosavljevic and Kotliar— formulated a variant of the 
DMFT where the geometrically averaged LDOS is com- 
puted from the solutions of the self-consistent stochastic 
DMFT equations. Subsequently, Dobrosavljevic et alm^ 
incorporated the geometrically averaged LDOS into the 
self-consistency cycle and thereby derived a mean-field 
theory of Anderson localization which reproduces many 
of the expected features of the disorder-driven MIT for 
non-interacting electrons. This scheme uses only one- 
particle quantities and is therefore easily incorporated 
into the DMFT for disordered electrons in the presence 
of phonons^ or Coulomb correlations^ In particular, 
the non-magnetic ground state phase diagram of the 
Anderson-Hubbard model at half-filling was derived.^ 

In this paper we investigate the Falicov-Kimball 
modcl^^ with a local disorder. The pure Falicov-Kimball 
model describes two species of particles, mobile and im- 
mobile, which interact with each other when both are 
on the same lattice site^SSi^ The Falicov-Kimball model 
captures some aspects of the Mott-Hubbard MIT, i.e. 
upon increasing the interaction the LDOS for mobile par- 
ticles splits into two subbands opening a correlation gap 
at the Fermi level if ng-l-n/ — 1, where rig (n/) is the den- 
sity of mobile (immobile) particles iSSiSL22iSS Here we in- 
troduce the Anderson-Falicov-Kimball model where the 
mobile particles are disturbed by a local random poten- 
tial. Our aim is to obtain a phase diagram of such a 
model and to identify MITs driven by correlations and 
disorder. We find a subtle competition between inter- 
action and disorder yielding stabilization of metalicity 
in the Anderson-Falicov-Kimball model. The model is 
solved within the DMFT framework combined with geo- 
metric averaging of the LDOS. The results are compared 
with those obtained within the DMFT with arithmetic 
averages. Only geometric averaging yields the Anderson 
transition. 
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In Section II we define the Anderson-Falicov-Kimball 
model and present the DMFT equations which provide 
the solution of this model. In Section III and IV numeri- 
cal results concerning the ground state phase diagram are 
shown and discussed in details. The analytical approach 
to the MIT in the band center is developed in Section V. 
Section VI presents our conclusions and final remarks. 



Bethe lattice with the band-width W the critical in- 
teraction obtained within the DMFT is Uc = W/2 for 
rie = n-f = i/2i^Si2LS2iSS How the disorder changes this 
gap and how the localized states enter into the mobile 
particle band are the subjects of the present study. We 
neglect any long-range order, which might be achieved 
on a fully frustrated lattice>i^ 



II. ANDERSON-FALICOV-KIMBALL MODEL 



B. Dynamical Mean— Field Theory 



A. The Model 

The Anderson-Falicov-Kimball model is defined by 
the following Hamiltonian 



V cj Cj 



E 



(1) 



where c\ (//) and Ci (fi) are fermionic creation and an- 
nihilation operators for mobile (immobile) particles at a 
lattice site i. Uj is a hoping amplitude for mobile parti- 
cles between sites i and j, and U is the local interaction 
energy between mobile and immobile particles occupying 
the same site. The ionic energy is a random, indepen- 
dent variable in our problem, describing the local dis- 
order disturbing a motion of mobile particles. We have 
assumed that only mobile particles are subjected to the 
structural disorder. This assumption could be relaxed in 
further generalizations of the model. 

The position of the immobile particles on a lattice is 
random if there is no long-range order. Therefore, we 
assume that the occupation number fjfi on the i-th site 
is equal to one with probability w {0 < w < 1) and zero 
with probability l — w. The presence of the randomly dis- 
tributed immobile particles introduces additional disor- 
der apart of that given by the e^-term in the Hamiltonian 
(^. However, the J7-term in the Hamiltonian |^ must 
be treated differently from the e^-term. The [/-term is 
operator valued in the immobile particle Fock subspace 
and one has to take the quantum mechanical average over 
a given quantum state of the /-particles. In contrast, the 
e^-term does not depend on the / operators and one has 
to average the quantum mechanical expectation values 
over different realizations of et or one has to study the 
whole statistics of an interesting quantum mechanical ex- 
pectation value. Whereas, the C/-term does not change 
the extended states into the localized ones, the e^-term 
can lead to such a change and thereby to the Anderson 
localization. 

In this paper we use the canonical description, where 
the number of the immobile and mobile particles are in- 
dependent of each other and fixed. In the pure Falicov- 
Kimball model increasing the interaction leads to open- 
ing a correlation (Mott) gap in the spectrum at i7 = 
(y^|26j27j28j29 If ^ = 1 the Fermi energy for mo- 
bile particles is inside of this correlation gap, which 
means that the ground state is incompressible. On the 



The Anderson-Falicov-Kimball model, where the in- 
teraction and disorder are local, is solved within the 
DMFT>i2iS242ii2£ Introducing a single and double-particle 
(Zubarev) Green functions^ Gij{uj) = ((ci|cj))„ and 

Tij{uj) — {{f^ fiCi\c^j))u,, respectively, we derive the equa- 
tions of motion 



j 



(2a) 



{u! + fj. - ek - U)Tki(uj) = {flfk)5ki + '^tkjT ji{uj) 
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(2b) 



where we used that a number of the immobile particles 
is conserved being zero or one, and hence flfifjfi — 
f-fi- The chemical potential /Lt is introduced only for 
the mobile subsystem and u> denotes the energy. Defin- 
ing the site-independent self-energy, according to DMFT 
scheme i^2i2S 



Gy [UJ) 



(3) 



which depends implicitly on e^, the system of the equa- 
tions Q can be solved yielding an explicite formula for 
the self-energy 



W(l - w)C/2 



Lu + fi ~ ei — {I — w)U — ri{uj) 



(4) 



We defined the averaged number of localized particles per 
site (f-fi) — Hf — w and we introduced the hybridization 
function ri{uj) which is a dynamical mean field [molec- 
ular field) describing the coupling of a selected lattice 
site with a rest of the system. The local non-interacting 
Green function is G°^{LL>) = l/[ti> + ^ — €i — r/(aj)] = 
GO(u;,e.). 

Using the self-energy E(a;, e^) and the hybridization 
function r]{uj) we obtain a local (e^-dependent) Green 
function 



Gii(w) = 



1 



+ /i - Ej - ri{uj) - T,{uj, Ej) 
and hence the ei-dependent LDOS 

1 



= G(cj,e,), (5) 



A(w,e,) 



-ImG(ti;, Ei). 



(6) 
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From the Ci-dependent LDOS Q we obtain either the 
geometrically averaged LDOS 

(w) = exp [(In ej))dis] (7) 

or the arithmetically averaged LDOS 

^arithlt^) = (^('^,e,:)>dis, (8) 

where (0(ei))dis — J deiP{ei)0{ei) is an arithmetic mean 
of O(ei).'^^ Here we used that are independent ran- 
dom variables characterized by a probability distribution 
function Vici). The lattice, i.e. translationally invari- 
ant, Green function is given by the corresponding Hilbert 
transform 

G(a.) = / (9) 

where the subscript a stands for either "geom" or "arith" . 
The Dyson self-energy S(a;) is determined from the k- 
integrated Dyson equation S(u;) = u) — rjioj) — 1/G(aj). 
The self-consistent DMFT equations are closed through 
the Hilbert transform G(w) = / deNo{e)/ [w - e - 
which relates the lattice Green function to the self- 
energy; here NQ{e) is the non-interacting density of 
states. 

The DMFT which uses Aarith is an exact approach in 
the limit of infinite dimension where quantum mechani- 
cal rescaling is imposed on hopping amplitudes^SiS^ On 
the other hand, the mathematically rigorous limit for the 
DMFT with Agcom is not yet known>i& Nevertheless, it 
is very promising single site theory which has the ability 
to describe at least some aspect of the Anderson local- 
ization, i.e. the localization due to random fluctuations 
of the wave function amplitudei24 As a single site theory, 
it cannot capture interference effects due to fluctuations 
of a phase of the wave function and, thereby, weak local- 
ization aspects are not recovered. 

The Anderson-Falicov-Kimball model Q is solved 
for a semi-elliptic density of st ates for the Bethe lat- 
tice, Noie) = iy^l- 4:(e/Wy/iTTW). Then r/(w) = 
W'^G{u!)/16. For a probability distribution function 
■p(ei) we assume a box model, i.e. V{ei) = 9(A/2 — 
|ej|)/A, with Q as the step function. The parameter A 
is a measure of the disorder strength. For numerical in- 
tegrations we use discrete values of Ci selected according 
to the Gauss-Legendre algorithm. The number of lev- 
els depends on A and is adjusted such to obtain smooth 
density of states. The chemical potential fj. ~ U /2, corre- 
sponding to half-filled conducting band (i.e., Ue = 1/2), 
and w = 1/2 are assumed in this paper. W = 1 sets the 
energy units. 

C. Criteria for Anderson and Mott MIT 

The arithmetically averaged LDOS Aarith(w) at the en- 
ergy u; in a band is always positive for non-interacting 
systems with disorder sii This quantity is non-critical for 



the Anderson localization. However it approaches zero 
when the gap is opened at energy to in the spectrum of the 
correlated system. On the other hand Agoom(^) vanishes 
at the Anderson localization. We therefore classify the 
states at energy lo as localized by disorder if Agoom('-^) = 
and Aarith('^) > 0. When both ylgeom(^) = and 
^arith('^) = it mcaus that the states at energy to are 
absent due to correlations. 



III. SPECTRAL PHASE DIAGRAMS 

Depending on the spectral properties of the Anderson- 
Falicov-Kimball model we distinguish three different 
regimes: i) weak interaction regime for < ?7 < W/2, ii) 
intermediate interaction regime for W/2 < U < 1.36W^, 
and iii) strong interaction regime for U > 1.36M^. Ex- 
amples of the spectral phase diagrams on the energy- 
disorder {u! — A) planes in these three regimes are shown 
in Fig. ^ in the upper, middle, and lower panels, respec- 
tively. 

In the weak interaction regime (i) the Mott gap is not 
opened. Increasing the disorder strength A leads to nar- 
rowing of the spectrum with extended gapless states (con- 
tinuous spectrum) and to broadening of the total band- 
width. This is illustrated in Fig. 13 presenting the evolu- 
tion of AgGom(w) and Aarith(w), Upper and lower panels 
respectively, upon increasing A. The continuous spec- 
trum corresponds to a support of ylgoom('j-'), i.e. such an 
energy window for which Agoom(w) > 0, whereas the full 
band is given by the support of Aarith(w). The trajecto- 
ries of the band edges (dashed lines) determined within 
the DMFT with arithmetic averaging and the trajec- 
tory of mobility edge (solid lines) determined within the 
DMFT with geometric averaging are shown in Fig.^ In 
the weak interaction regime the localized gapless states 
(pure point spectrum) are in a compact part of the spec- 
tral phase diagram between the mobility and band edges. 
The spectral phase diagram of the Anderson-Falicov- 
Kimball model in the weak interaction regime (i) is qual- 
itatively similar to that of the Anderson model without 
the interactionJ:^: 

In the intermediate interaction regime (ii) the Mott 
gap is opened at A = 0, as is shown in Fig. El Upon 
increasing A the Mott gap is shrunk and finally closed, 
cf. Fig. 13 The disorder redistribute the spectral weight 
filling in the correlation gap completely by continuous 
spectrum. The spectral phase diagram is shown in the 
middle panel of Fig. As previously the total band- 
width, determined by ylaiith('^), increases by increasing 
A. However in contrast to (i), in the present case there 
are two trajectories representing external (from the band 
gap sides) and internal (from the Mott gap sides) mo- 
bility edges. Similarly there are two band-edge trajec- 
tories, external and internal ones. The extended gapless 
states are bounded between mobility edge trajectories as 
is shown in the middle panel of Fig. ^ Two separated 
regions with localized gapless states are bounded by the 
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mobility edge and the band edge trajectories. If the cor- 
relation gap is opened, the localized states appear in the 
spectrum from each side of the Hubbard subbands. 

In the strong interaction regime (iii) the Mott gap, de- 
termined within the DMFT framework with geometrical 
averaging, is never filled in even at large disorder. The 
LDOS given by Agcom(w) and shown in Fig. 0] always 
has two separate parts, remnants of the lower and the 
upper Hubbard subbands. This is in contrast with the 
DMFT framework with arithmetic averaging where these 
two subbands always merge if A is sufficiently large. It 
means that the Mott gap is only filled in by localized 
states when the disorder increases. In the spectral phase 
diagram the spectrum of extended states is given by two 
separate lobes bounded by the mobility edge trajectories, 
as in the lower panel of Fig. ^ At A > these lobes are 
surrounded by localized gapless states. 



IV. BAND CENTER 

In the half-filled band case the ground state proper- 
ties are solely determined by the character of quantum 
states in the band center {uj = 0). Corresponding phase 
diagrams in the interaction-disorder {U — A) plane are 
shown in Fig. [S] and discussed below. We define three 
phases: a) extended gapless phase (i.e. a gapless phase 
with extended states at the Fermi level), b) localized gap- 
less phase (i.e. a gapless phase with localized states at 
the Fermi level), and c) gapped phase (i.e. a phase with 
a gap at the Fermi level). 

The extended gapless phase (disorder metallic phase), 
is characterized by a non-zero value of the LDOS. In 
the pure Falicov-Kimball model the Luttinger theorem 
is not satisfied and quasiparticles at the Fermi level have 
finite life-time. ■^^ It means that due to the interaction 
Aa{uj = 0) < No{u! ~ 0) even at A = 0. Increasing A 
at fixed U leads to further decreasing of the LDOS as is 
shown in the upper panel of Fig. Similarly, increasing 
U at constant A leads to decreasing of the LDOS as is 
presented in the upper panel of Fig. 

Mott-Hubbard MIT, represented by Af"{U) lines in 
Fig. [SJ occurs at small and intermediate disorder < 
A < 1.70W and the interaction W/2 < U < 1.36W. 
This MIT is continuous one as is seen in the lower panel 
of Fig. El around A « 0.8 and in Fig. □ around U « 0.7 
(upper panel) and U ~ 1.3 (lower panel). Increasing 
A above Af ^(f7) when W/2 < U < 1.36W results in 
a transition from a correlated gapped insulator into a 
bad metallic phase, as is illustrated in the lower panel of 
Fig. It means that the disorder stabilizes the metallic 
phase. Similar Mott-Hubbard MIT is obtained within 
the DMFT with arithmetic averaging and the corre- 
sponding phase diagram is reproduced in the lower panel 
of Fig.El 

Anderson MIT fine A^{U) is an increasing function 
of U for < U < 0.95W starting from the value 
A^ = eW/2 « 1.361^ (e « 2.718 is the Euler constant) in 



the non-interacting casein The interaction impedes the 
localization of particles due to impurity scattering. In 
particular, for eW/2 < A < 2mW and < t/ < 0.951^ 
the interaction turns the Anderson insulator into a bad 
correlated metal, as is shown in the lower panel of Fig.|7| 

Mott and Anderson insulators are rigorously defined 
only for U > W/2 with A = and only for A > eW/2 
with [/ = 0, respectively. In the presence of the interac- 
tion and disorder this distinction can no longer be made. 
However, as long as the LDOS shows the characteristic 
Hubbard subbands (see Fig. ^ one may refer to a dis- 
ordered Mott insulator (gapped phase). With increasing 
A the spectral weight of the Hubbard subbands vanishes 
and the system becomes a correlated Anderson insula- 
tor (localized gapless phase). The border between these 
two types of insulators occurs at A{U) « eW/2V2 when 
U ^ W. To estimate this value we used the analyti- 
cal result from Ref. ild Our estimation is exact when 
U ^ oo since the band-width of the Hubbard satellites 
is rigorously known to be W/\f2 for the Falicov-Kimball 
modcL^ The results obtained here within DMFT show 
that the Mott and Anderson insulators are continuously 
connected. Hence, by changing U and A it is possible to 
move from one type of the insulator to the other with- 
out crossing the metallic phase. This is plausible because 
the Anderson MIT {U = 0) and the Mott-Hubbard MIT 
(A = 0) are not associated with a symmetry breaking. 



V. LINEARIZED DYNAMICAL MEAN-FIELD 
THEORY 



At the MIT (dots and squares in Fig. ISJ the LDOS 
vanishes in the band center. Therefore, in the vicinity 
of the MIT but on the metallic side the LDOS is arbi- 
trary small and the transition points on the phase dia- 
gram can be determined analytically by linearizing the 
DMFT equations >iSi2i In the band center, due to a sym- 
metry of Ac((a;), we find that 0(0) — —it^AaifS) and 
is purely imaginary. Hence the DMFT self-consistency 
leads to the following recursive relation ?7("+i)(0) = 

-inW'^J^a^ (0)/16, where the left hand side in the (n-|-l)- 
th iteration step is given by the result from (n)-th itera- 
tion step. Using Eqs.|Sl with ^ and expanding them with 
respect to small ^i"''(0) we find from Eq. Elthat 



^("+i)(0,6,) = ^A(r)(0)T(60, 



(10) 



where 



T(6) 



-(If 



(11) 
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The recursive relations within the linearized DMFT (L- 
DMFT) with geometrical averaging are 



^(«+i)(0) = (0) 

geom V ^ / geom V ^ / 



1 f^'^ 

- / delnT(6) 

^ J -A/2 



(12) 



and within L-DMFT with arithmetical averaging are 



arith (0) = ^arithC^)^^ 



A/2 



rfeT(e) 



-A/2 



(13) 



In a metallic phase the recursions are increasing, i.e. 
aL"+'^(0) > Ai"H0), wher eas in the insulating phase 
they are decreasing. Therefore, at the boundary curves 
between metallic and insulating solutions in Fig.|Slthe re- 
cursions are constant Aq""'^^'' (0) = Ai"''(0). This observa- 
tion leads directly to the exact (within DMFT) equations 
determining the curves A = A(J7), i.e. 



1 = exp 

16 ^ 



A/2 



-A/2 



delnT(e) 



exp [/g, 



for L-DMFT with geometrical averaging, and 



.(C/,A)] 
(14) 



16 



A/2 



deT(e) 



-A/2 



16 



/anth(f/,A), (15) 



for L-DMFT with arithmetical averaging. Both integrals 
can be evaluated analytically with the results 



In 









t) 


- 2 In 




2U 








arctan 











;com(t/,A) = 
2 / A \ 2' 



A 
2" 

A- 



U 



U 



and 



rith 



(f/,A) 



(16) 



(17) 



Solutions of Eas. 1141 and 1151 are shown as solid curves in 
the upper and the lower panels of Fig. respectively. 
We find excellent agreement with the numerical solu- 
tions of the full DMFT equations (dots and squares in 
Fig. 13). At small U the critical disorder strength ob- 
tained from (|14|l increases linearly with the interaction, 
i.e. A(C/) « We/2 + ttU/2. This is because the total 
bandwidth increases linearly with U . At small A the so- 
lution of Eq. ^] (L-DMFT with geometrical averaging) 
is found to be A([/) w ./iP 
the curve A(C/) 



{W/2Y. It turns out that 



yJU"^ — {W/2Y is also an exact solu- 
tion of Eq. El (L-DMFT with arithmetical averaging) for 



all U > W/2. At weak disorder both approaches give the 
same results. 

VI. CONCLUSIONS 

In the present paper we introduced the Anderson- 
Falicov-Kimball model and solved it obtaining the phase 
diagrams for the ground state with the suppressed long- 
range order. Three different phases, the disordered 
metal, the disordered Mott insulator, and the correlated 
Anderson insulator, were identified. It was shown that 
correlation and disorder compete with each other sta- 
bilizing the metallic phase against occurring one of the 
insulators. We found that these two insulators are conti- 
nously connected. 

The phase diagram with the three phases in the ground 
state for the Anderson-Falicov-Kimball model is similar 
to the phase diagram for the Anderson-Hubbard model 
solved within the DMFT with geometric averaging in 
Ref.H There are however important qualitative differ- 
ences between these two solutions. In the Anderson- 
Falicov-Kimball model the Mott transition is continuous 
whereas in the Anderson-Hubbard model there is a hys- 
teresis at low and a crossover transition at high disorder 
strengths. In addition, in the Anderson-Hubbard model 
the Luttinger pinning in the disordered metal is recon- 
structed by strong correlations. This feature is absent 
in the Anderson-Falicov-Kimball model, where the Lut- 
tinger pinning is violated even in the pure case*^ 

A similar technique, i.e. the DMFT with geometric av- 
eraging, could be used to solve other models with disor- 
der and interaction between quantum (mobile) and clas- 
sical (immobile) degrees of freedomi^ In such cases the 
self-energy should be given analytically and this removes 
the problem of using any numerical impurity solver, as 
was necessary in the Anderson-Hubbard model. ^ Then 
not only the LDOS at the Fermi level but also mo- 
bility and band edge trajectories can be easily deter- 
mined. Such models are important for understand- 
ing the physics of manganites'^^ or diluted magnetic 
semiconductors where charge carriers are coupled to 
randomly distributed localized magnetic moments. The 
role of the disorder and the Anderson localization are 
inherent for those correlated systems. 
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FIG. 1: Spectral phase diagrams of the ground states without 
a long-range orders for the Andcrsori-Falicov-Kimball model 
with w — 0.5 at f/ = 0.3 (upper panel), U — 0.9 (middle 
panel), and U = 1.5 (lower panel). Solid lines show mobility 
edges determined within DMFT with geometric averaging and 
dashed lines present band edges determined within DMFT 
with arithmetic averaging. 
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FIG. 2: Geometrically (upper panel) and arithmetically 
(lower panel) averaged local density of states at w = 0.5 and 
U = 0.3 for different disorder strength A. Vanishing of Ageom 
upon increasing A is a signature of the Anderson localization. 




FIG. 3: The same as in Fig. El but at [/ = 0.9. The correla- 
tion (Mott) gap is closed by strong disorder and the system 
becomes a bad metal. 




FIG. 4: Geometrically averaged local density of states at w = 
0.5 and U = 1.5 for different disorder strength A. 




FIG. 5: Ground state phase diagrams for particles in a band 
center determined by using geometric (upper panel) and axith- 
metic (lower panel) means. Dots and squares are determined 

from the numerical solution of the DMFT equations. Solid 
lines are obtained analytically from the linearized DMFT. 




FIG. 6: Local density of states in a band center (cj = 0) as 
a function of disorder A at to = 0.5 with U = 0.3 (upper 
panel) and U = 0.9 (lower panel). Solid (dashed) lines are 
determined by using geometric (arithmetic) averaging. 




FIG. 7: Local density of states in a band center (cj = 0) as 
a function of interaction U a.t w = 0.5 with A = 0.6 (upper 
panel) and A = 1.6 (lower panel). Solid (dashed) lines are 
determined by using geometric (arithmetic) averaging. 



